##Replication code for "drivers of resemblance: explaining the conversion of pre-electoral coalitions
##into coalition cabinets"
library(QCA)
library(tidyverse)
library(ragg)
library(SetMethods)
library(showtext)

sec_ch <- read.csv("dpr.csv", sep = ";") %>%
          filter(type == "pre") %>%
          select(-type, -legislative_type, -inter_party_primaries,
                 -low_chamber, -upp_chamber, -status, -ce) %>%
          filter(country != "Chile")

sec_ch1 <- sec_ch %>%
           column_to_rownames(var = "cases") %>%
           select(-country)

sec_ch2 <- sec_ch1 %>%
           mutate(lwpol = calibrate(lwpol, type = "fuzzy", method = "direct", threshold = "e=3, c=2, i=1")) %>%
           mutate(hlpol = calibrate(hlpol, type = "fuzzy", method = "direct", threshold = "e=1.5, c=2.7, i=4")) %>%
           mutate(htemp = calibrate(htemp, type = "fuzzy", method = "direct", threshold = "e=85, c=70, i=55")) %>%
           mutate(lplp = calibrate(lplp, type = "fuzzy", method = "direct", threshold = "e=0.5, c=0.405, i=0.3")) %>%
           rename(LWPOL = lwpol,
                  MAJ = maj, 
                  HLPOL = hlpol,
                  HTEMP = htemp,
                  LPP = lplp,
                  CR = cr)

##Analysis of Sufficiency
TTsec1 <- truthTable(sec_ch2, outcome = "CR", conditions = c("MAJ", "LWPOL", "HLPOL", "HTEMP", "LPP"),
                     complete = TRUE, incl.cut = 0.8,
                     sort.by = "incl, n+",
                     show.cases = TRUE)

ESA <- findRows("~MAJ*~HTEMP", obj = TTsec1)

##The code below produces Table A.8 in the Supplementary Material. 
intermediate_result <- minimize(TTsec1, include = "?", details = TRUE, dir.exp = "1,1,1,1,1", exclude = ESA)

stargazerSol(intermediate_result, 
             outcome = "CR",
             show.cases = T,
             type = "latex")

